
SUBSTITUTE SPECIFICATION 
METHOD FOR THE LINEARIZATION OF FMCW RADAR DEVICES 



CROSS REFERENCE TO RELATED APPLICATIONS 
This application is based upon German Application 
Serial No. 10315012.9 Filed on April 2, 2003 wherein 
priority is claimed under 35 U.S.C. 119. 

BACKGROUND OF THE INVENTION 

The invention relates to a method for the linearization 
of a frequency modulated continuous wave (FMCW) radar 
device . 

With FMCW radar devices, phase errors occur in the 
reception signal, as the result of a not strictly linear 
characteristic in a transmission frequency of an FMCW radar 
device. This results in a loss of accuracy in the 
evaluation of the signal. 

In this case, various methods are known for linearizing 
an FMCW radar device. For example, one known method 
discloses a method wherein a delayed transmission signal is 
generated via an optical delay line having a defined running 
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time. In this case, the delayed transmission signal is 
mixed with the non-delayed transmission signal. Thus, it is 
possible to determine the deviation from the ideal linear 
characteristics. In particular, the disadvantages of this 
method are the great amount of construction effort, expense 
and the age-related clouding of the optical components. 

In addition, it is also known to generate a signal 
function that can be predetermined, such as a linear 
characteristic via a digital direct synthesis (DDS) system. 
However, a disadvantage of this system is that there is a 
large construction effort and expense associated with this 
cost intensive structure. 

An FMCW radar system having linear frequency modulation 
is known from German Patent DE 4104792 Al . In this system, 
there is an oscillator that generates a ramp- shaped 
modulated signal having a modulated oscillation frequency, 
which is emitted by an antenna. A monitoring and control 
unit determines the deviation of the measured modulated 
oscillator frequency of a signal reflected from an object 
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via a Hilbert transformation using a linearized reference 
value. In this case, a phase deviation can be determined 
from the deviation that has been determined. 

From this deviation, a correction signal is generated, 
which is passed to the oscillator, whereby the oscillation 
frequency is modified in accordance with the predetermined 
reference values. A disadvantage of this system is that 
this involves a great construction effort, which is cost 
intensive, since additional components are required. 

SUMMARY OF THE INVENTION 

The invention relates to a method with which an FMCW 
radar device can be linearized, without a great construction 
effort or expense. 

In this case, a correction phase term is calculated on 
the receiver side for compensation of the phase error in the 
reception signal q(t) . In contrast with known methods, no 
costly components are required in the transmission branch of 
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the FMCW radar device when using this method and this thus 
results in reducing costs. 

This method, which relates to the invention should 
include the following steps: 

selecting a number L of consecutive ramp- shaped 
reception sequences q k (n) of the reception signal, which 
number can be predetermined, with K=1,...,L, the phases 
arg{q k (n)} of which can be represented as a polynomial of the 
N th order for the time index n, with the polynomial 
coefficient m*, with {!=!,..., N, 

transforming a spectrum range Q(e jQ )of the selected 
rece ption sequences q(n) into a basic band reception 
sequence q k {n) with k=0,...L-l, 

iteratively calculating a correction phase term for 
partial compensation of non-linear frequency components in 
the basic band reception sequences q k (n) by means of 
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calculating polynomial coefficients of the individual 
basic band reception sequences q k (n) via estimating methods 

wherein ^( n ) are the sequences that have already been 
iteratively phase-corrected, and whereby the iteration is 
stopped once the parameter change between two consecutive 
iterations, which can be predetermined, remains below a 
threshold e , which can be predetermined. 

In another embodiment of the invention, the scanning 
cycle T A of the ramp signals is reduced, in a method step 
after the base band transformation. For this purpose, a 
prior antialiasing low pass filtering of the ramp signals 
q k (n) is performed. In addition, it is beneficial if the 
downsampling factor K, by which the scanning cycle T A is 
reduced, lies between K=30 and K= 60. 

In contrast, with an alternative linearization of the 
frequency modulation of the radar transmission signal, the 
linearization is performed on the receiver side, via a 
digital processing of the reception signal. Thus, no 
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additional components are required for this method. The 
calculations required for this method can be performed by 
the components for signal evaluation that are already 
present. This results in additional advantages with regard 
to the production costs. 

Another advantage of the invention is that there can be 
in an inexpensive estimation method that can be easily 
implemented, since the reception signal is processed in the 
base band. This signal results in a reduction of the 
scanning cycle T A by a factor K between 30 and 60. In 
addition another advantage of this method is that there can 
be a calculation of the polynomial coefficients via a robust 
estimation method. This is even in the case of channels 
subject to severe interference. Thus, there can be a 
relatively good compensation of phase errors with signal to 
noise ratios (SNR) of up to -20 dB using this method. Thus, 
the method can be used for not only the civil sector, but 
also in the military sector or for distance measurements. 

BRIEF DESCRIPTION OF THE DRAWINGS 
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Other objects and features of the present invention 
will become apparent from the following detailed description 
considered in connection with the accompanying drawings, 
which discloses at least one embodiment of the present 
invention. It should be understood, however, that the 
drawings are designed for the purpose of illustration only 
and not as a definition of the limits of the invention. 

In the drawings, wherein similar reference characters 
denote similar elements throughout the several views: 

FIG. 1 is an exemplary idealized ramp-shaped FMCW 
transmission signal; 

FIG. 2 is an exemplary representation of the refection 
of the radar signal on a point shaped object at distance R; 

FIG. 3 is a schematic diagram of the method that 
relates to the invention; 
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FIG. 4 is a representation of the base band shift of a 
selected spectrum range in the reception signal; 

FIG. 5a is a representation of a spectrum of a 
downsampled basic band signal before and after the band pass 
filtering; and 

FIG. 5b is a representation of a magnitude spectrum of 
a radar reception signal before and after the phase 
correction . 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Turning now in detail to the drawings, FIG. 1 shows the 
starting point of this method is the transmission signal of 
an FMCW radar that is in the general polynomial approach of 
the N th order for a transmission ramp as: 

x{t) ^ xvU 2r{fj + ^ + ^ + ^ + ... + r 0 )} (1) 

With /e[0,r r ]. With this design, the coefficients m 
with 1-2,... ,N represent the polynomial parameters, yo 
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represents the initial phase, and T r represents the time 
duration of the transmission ramp. The parameter m x or slope 
that is assumed to be given, indicates the linear component 
mi=B/t # which is exclusively present in an ideal case, 
wherein B describes the frequency deviation of the FMCW 
radar system. Thus in the following formula, a single 
reflecting point shaped object is assumed. Thus, the 
corresponding parameters are represented in FIG. 1, using 
the example of an idealized transmission signal. Thus with 
Equation (1) the reception signal in an ideal, noise free 
case is then: 

y{t) = tx V {jlr\f 0 t^-t 2 +-t> +_r 4 +...-f /o -—J) (2) 

Since the determination of relative distances between 
objects that are moving relative to the radar sensor at 
approximately the same speed represents a possible use of 
the method described, it can be assumed that s(t) = 2R in 
the following (Fig. 2) . 
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For the momentary wavelength Mt) , the following 
results at the speed of light c: 



l- = £& = ±L+±mA with f(0 = 

(t) C C\ m J 



1 darg{x(t)} 
In dt 



/o+I>/ 



N 



It is advantageous if the reception signal y(t) is mixed 
into a lower frequency position. The mixed reception signal 
q e (t) results as follows, from Equations (1) and (2): 



If q e (t) describes a harmonic exponential function, i.e. 
if m f s 0 with 0 = 2, . . . , N, then the ideal case is present, 
and the transmission signal x(t) possesses a purely linear 
increase in the transmission frequency. 

In order to indicate the maximal deviation from this 
ideal transmission frequency progression, a maximal relative 
linear error is defined as 



9,(0=^(0 'X (0= exp \- j—— [fo + L™** 




(3) 
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1/(0 -f ideal (Q\ 

0:= max — 



i 

= — max 
/o-0 w, 



£=2 



with f id (t)=f 0 +mit for the ideal case. The symbol 0 
therefore describes a maximum deviation of the momentary 
frequency progression for f 0 = 0 from the linear function 

f..(t) I . Taking into consideration the noise term r(t) that 

J id V / l/ 0 =0 

is always present in real applications, the following is 
obtained from Equation 3, with q(t) = q e (t)+r(t) 



g(0= expj- J^f\fo + f + KO (4) 



If the complex-value noise signal r(t) is assumed to be 
white and gaussian distributed in Equation 4, for the sake 
of simplicity, the noise influence for a significantly large 
signal-to-noise ratio can also be considered as phase noise 
e j7{,) , whereby f{t) is real-valued, white, and gaussian 
distributed. The phase signal, which is made time-discrete 
with the scanning cycle T A , results from 
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arg 



{q(nT A )}=4(nT A ) 



c 




(5) 



Proceeding from the known (received) phase scanning 
values <j)(nT A ), the method according to the invention 
implements a robust estimate of the polynomial parameters mj, 
with t = 2, . .., N. In an advantageous embodiment of the 
invention, a robust estimate of the unknown object distance 
R is additionally possible. 

Fig. 3 shows a block schematic in the form of a 
sequence schematic of the method according to the invention. 
The initialization 1 takes place, according to the 
invention, by means of a base band shift BV of L 
consecutive ramps of the sequence q(n) = q(nT A ) that results 

from Equation (4) , by means of time discretitation . The 
shift of a spectrum range, which can be predetermined, of 



Q(e jCl> ) = F {<i{ n )} into the base band is shown as an example in 

Fig. 4. The cut-off frequencies Qi and Q u are defined, for 
example, by way of the minimum and maximum distances R m i n and 





12 



Rmax/ which are determined by the purpose of use of the FMCW 
radar system: 

A <TT A Tf j 

— R^rn,T A , Cl u =—R^m x T A , Q,„=-(Q„-Q,) 
c c l 

The ramp signals q k (n) that are formed in this manner 
possess a band-pass spectrum. It is advantageous that after 
the base band shift and the subsequent antialiasing low- 
pass filtering TF, the scanning cycle can be reduced by a 
factor K. In this connection, the factor K particularly 
lies between K = 30 and K = 60. This results in advantages 
with regard to a reduction in effort/expense of all further 
partial operations . 

The next step or progression of the block diagram in 
Fig. 3 shows the sequence of the iterations for the 
calculation of a correction phase term for partial 
compensation of non-linear frequency components in the basic 
band reception sequences q { p(m) with m as a time index after 

cycle reduction. In the further process steps, a signal 
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sequence whose scanning cycle T A was reduced is assumed with 
qf(m) - Of course, these process steps can also be carried 
out with a sequence qf(n), whose scanning cycle T A was not 
reduced. 

Within the loop that runs through in each iteration, it 
is advantageous if band-pass filtering of the sequences 

qf(m) is located at the beginning, which is advantageously 
carried out as windowing in the spectrum range, whereby all 
the spectral coefficients outside the window are set to 
zero. It is advantageous if a Hamming window is used as the 
window function. 

In the case of several reflecting objects, the 
reception signal spectrum of a single object is selected by 
means of the band-pass filtering. All of the other, 
undesired spectra are suppressed. In addition, suppression 
of those noise -like spectrum components that do not belong 
to the working component takes place. In this way, the SNR 
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can be increased after the band-pass filtering, even in the 
case of a single reflecting object. 

To calculate the band-pass filtering, it is 
advantageous to perform a Fast Fourier Transformation FFT 
of the sequences qf{m) . For this purpose, it is 
advantageous if the individual discrete Fourier transforms 
Ql'\ju) of the base band reception sequences g[ n (m) are 
calculated, whereby Ql n (M) = FFT{q ( ?(mj\ for k = 1,..., L. The 

number of FFT points M typically lies at up to 1000. In the 
case of a reduction of the scanning cycle T A that has 
advantageously been performed, and prior antialiasing low- 
pass filtering, the FFT can be carried out with a smaller 
number of points M of 64 to 256. 

In a next iteration step, filtered base band reception 
sequences q ( ^{m) are calculated by means of a band-pass 

filter according to Q?\ju)= w(/j) Q[ n (ju) , with qU) ^ = FFT{q ( k ° (m)} 
, whereby w(u) is a spectrum window that can be 
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predetermined, and //€[//,, indicates the range of the 

spectrum window having a window center point i^ax that can be 

predetermined, whereby, with a lower limit \i u that can be 

predetermined, and an upper limit |a i that can be 
predetermined. In order to determine the position of the 
center point |i ma x, the FFT magnitude spectra Q^\jJ) 



are 



averaged over the number L, and the maximum is sought from 
the averaged spectrum, according to the following formula: 



G (0 (/O = max -Y \q^(/i)\ 



For example, for a single reflecting object, \x i = 0 and 
ja u = M-l. For several non- over lapping objects, the limits (i i 
and \x u must be selected so that in each instance, only the 
spectral component that corresponds to the object of 
interest is detected. 
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Subsequently, it is advantageous if the calculation of 
the individual inverse discrete Fourier transformed q ( k l) (m) of 

the filtered FFT spectra Q k {i) (ju) takes place, whereby 
qj t i \m) = IFFT{Q^(ju)} for k = 1, . . . , L. 

In a subsequent iteration step, the distances Rk (l) are 
estimated by means of a maximum likelihood estimating 
method. The polynomial coefficients Rk 0) can be estimated 

from the estimated distances ™a . Subsequently, the 
distances R k are averaged over L reception sequences q k (n) ■ 
In a final iteration step, a reception sequence £< ,+,) (m) with 

the averaged estimated polynomial coefficients m i is 
calculated as the starting point for a next iteration loop, 
according to 

f d7rR U) n 1 
q ( r\m) = a°(«)-«p|./^-5j«l < (iii«rj l | 
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The maximum likelihood parameter estimate takes place by 
using the linear estimator 



v = [A T Ay A T b (6) 

with the vectors y = [Rf Q9 Rm } ,Rm 2 ,...] T , b = [Q(0)MT A )M2T A ),..] r and 

a matrix A having the dimension L q x (N+l) , whereby L q 
represents the length of a sequence q { p{m) . The 

determination of the phase values <t>(nT A ) in the vector b 
takes place by means of a phase unwrapping method. 

A reduced linearity error in the downsampled base 
band signal sequence q[ n (m) is successively achieved by means 

of the phase correction, and this results in a lesser 

estimate variance, in each instance. If the condition 

IP'" 1 * - j? (0 |< s has been fulfilled after i = I iterations, final 

estimated values R and ™e can be calculated from the 
individual intermediate results R K9 and , which have been 
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averaged over all L reception sequences in the i iteration, 
in each instance, according to: 

R = R U) and m t = ^£F>m<<> , /-2....N 
R 1=1 

In the following, it will be shown that it is possible, 
using the method according to the invention, to reconstruct 
the spectral resolution, which has been reduced as a result 
of the linearity error in the FMCW radar transmission 
signal, of the sampled, mixed reception signal q(n) . This 
will be shown in exemplary manner below, using the example 
of two objects at different distances. 

The starting point is a quadratic transmission signal 
sequence (chirp) with a linearity error of 0=5%, whereby 
the simulation parameters are selected to be N = 2, m 2 = 
5xl0 12 , SNR = -18 dB, K = 40, and M = 64 . For the other 
system parameters, the following table applies: 
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Parameter 


Value 


Minimum possible distance 


Rmin = 900 m 


Maximum possible distance 


Rmax = 1,100 III 


Ramp increase time 


T r = 1 ms 


Radio Fundamental frequency 


f 0 = 35 GHz 


Frequency deviation 


B = 100 MHz 


Sampling frequency 


f A = 2 MHz 


Spectrum averaging L 


L = 20 


Spectrum averaging, total 


25 x L = 500 


Spectrum window w(|a) 


Hamming window 



The objects are located at the distances R 0 = 900 m and R x = 



1,100 m. Fig. 5 shows the magnitude spectra 



\QM J %\Q n {e jn ) 



of the downsampled signal q ( k °\m) for a fixed k = 0, L-l 

during the first iteration before (top representation in 
Fig. 5a) and after (bottom representation in Fig. 5b) band- 
pass filtering. The partial spectra, which can be easily 
separated, are clearly evident, whereby the bottom partial 
spectrum, which belongs to the object at the distance R lt is 
selected in the band-pass filtering. 



The top figure in Fig. 5b shows the range of the 
agnitude spectrum \Q{ eJn ) of the radar reception signal that is 
of interest, which is established by the minimum and maximum 



m 
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distance (see table) . As the result of the non-linearity in 
the transmission signal, the frequency lines that correspond 
to the two objects are smeared and can hardly be recognized. 
Using the parameters estimated for the object at the 
distance Ri, using the method according to the invention, the 
inherently non- linear component of this signal is now 
compensated, as shown in the bottom part of Fig. 5b. 

A comparison of the figures shows the clear increase in the 
object resolution. The small triangles in the bottom part 
of Fig. 5b characterize the positions of the frequency lines 

in the mangitude spectrum \Q{ ejQ ) , which would generate a 
transmission signal having a linear frequency increase. 

Accordingly, while at least one embodiment of the 
present invention have been shown and described, it is to be 
understood that many changes and modifications may be made 
thereunto without departing from the spirit and scope of the 
invention as defined in the appended claims. 
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ABSTRACT OF THE DISCLOSURE 



A method for the linearization of frequency modulated 
continuous wave (FMCW) radar devices having non- linear, 
ramp-shaped, modulated transmitter frequency progression 
x(t) . With this invention, a correction phase term for 
compensation of the phase error in the reception signal q(t) 
is calculated on the receiver side in this device. 
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